The lateral Casimir force beyond the proximity force approximation : a nontrivial 
interplay between geometry and quantum vacuum 
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The lateral Casimir force between two corrugated metallic plates makes possible a study of the 
nontrivial interplay of geometry and Casimir effect appearing beyond the regime of validity of the 
Proximity Force Approximation (PFA). Quantitative evaluations can be obtained by using scattering 
theory in a perturbative expansion valid when the corrugation amplitudes are smaller than the three 
other length scales: the mean separation distance L of the plates, the corrugation period Ac and the 
plasma wavelength Ap. Within this perturbative expansion, evaluations are obtained for arbitrary 
relative values of L, Ac and Ap while limiting cases, some of them already known, are recovered 
when these values obey some specific orderings. The consequence of these results for comparison 
with existing experiments is discussed in the end of the paper. 



I. INTRODUCTION 

The Casimir effect [l| is the dominant interaction be- 
tween neutral plates separated by distances in the mi- 
cron or submicron range. The better and better con- 
trol of this Casimir effect achieved over the last 10 years 
0, H, 0, H 0, 0] hence opens new roads for the design of 
nanoelectromechanical systems (NEMS) Q. 

Casimir initially studied the simplest geometric con- 
figuration with two parallel plane plates large enough so 
that the theoretical analysis is simplified thanks to the 
lateral translation symmetry 0]. Except for a few cases 
0, experiments are performed between a plane and a 
sphere H, @, a geometry more easily mastered at dis- 
tances in the micron or submicron range. Force evalua- 
tions in this geometry are commonly obtained by using 
the so-called proximity- force approximation (PFA) [l(|, 
with the energy simply obtained by averaging the plane- 
plane expression over the distribution of local separation 
distances met in the plane-sphere geometry. It is com- 
monly agreed that this approximation is valid when the 
radius R of the sphere is much larger than the closest 
separation distance L 11(. But a quantitative determi- 



nation of its accuracy in plane-sphere experiments is still 
lacking, at least for the problem of experimental interest 
where electromagnetic fields are reflected on metals, the 
optical properties of which must be accounted for at the 
distances met in the experiments. In contrast, results 
valid beyond the PFA have been reported for theoretical 
models involving scalar fields reflected on perfect bound- 
ary conditions [TH, [IH, [l4| • 

At this point, we may stress that the situations which 
may be treated within the PFA correspond to a trivial 
interplay between geometry and Casimir effect since the 
geometry is described by an averaging over the distribu- 
tion of local distances. In contrast, the general case opens 
a far richer physics with a variety of stimulating theoret- 
ical predictions [HI, [H, [uj, so that the exploration of 
situations beyond reach of the PFA raises great expecta- 



tions. The idea can already be tackled for the description 
of the effect of roughness on the Casimir force. This de- 
scription is commonly given within the PFA [181 ] valid 
only when the wavelengths associated with the plate de- 
formation are large enough [l£| l20|, |2lJ. As the effect 
of roughness is only a small correction of the Casimir 
force, and the characterization of the roughness state of 
the plates is not very accurate, one can hardly expect 
quantitative theory-experiment comparisons in this case. 

Fortunately, there exists a geometry better suited to 
the aim of an accurate theory-experiment comparison, 
namely that with parallel and periodic corrugations im- 
printed on the metallic surfaces. The Casimir force con- 
tains a lateral component besides the usual normal one, 
since lateral translation symmetry is broken here [22| . 
The lateral Casimir force is smaller than the normal one, 
but it has nevertheless already been measured in exper- 
iments [23j]. It is easily computed within the PFA and 
has also been calculated beyond the PFA by using more 
elaborate theoretical methods. The lateral force has first 
been evaluated for perfect mirrors using a path-integral 
formulation in a perturbative [24j or non perturbative 
approach [2|| . As expected, the PFA is found to be valid 
only in the limiting case where the corrugated surfaces 
are nearly plane for the vacuum fields involved in the 
calculation of the Casimir energy. When introducing the 
corrugation wavelength Ac and the mean separation dis- 
tance L, ones characterizes this limit as Ac 3> L. 

Now, the experiments have been performed with dis- 
tances L around 200nm at which it is essential to take 
the optical properties of the metals into account [2(| [13] • 
A simple description of these optical properties is given 
by the plasma model introducing a further length scale, 
the plasma wavelength Ap, with a typical value of 137nm 
for Gold plates. A novel theoretical method has recently 
been presented [H| which allows one to calculate the lat- 
eral Casimir force for arbitrary relative values of the three 
length scales L, Ac and Ap, provided that the corrugation 
amplitudes a\ and ai are smaller. The approach relies 
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on scattering theory Q used in a perturbative expansion 
pl| with respect to the corrugations. 

In the present paper, we will first present the complete 
derivation of the results presented in (28[ . We will write 
the lateral Casimir force at the order oc aiai for arbi- 
trary relative values of L, Ac and Ap. Limiting cases will 
be obtained when these values obey some specific order- 
ings. Some of them are already known, in particular the 
PFA limit and the perfect reflection limit, and the known 
results are recovered as expected in our calculations. 

We will conclude the paper by discussing the con- 
sequence of these results for comparison with existing 
measurements [23]. These measurements were found to 
agree with the PFA computations, within the margins 
of experimental uncertainty [23[ . As this conclusion dif- 
fers from that drawn from our results [28j |. there may 
be a potential concern for theory-experiment compari- 
son, and the question will be discussed in detail at the 
end of this paper. Let us recall that the corrugation am- 
plitudes used in the experiments [23| were smaller, but 
not much smaller, than the other length scales, so that 
the theoretical predictions drawn from our perturbative 
expansion cannot be compared directly with the experi- 
mental results [29l [30j. However, it seems unlikely that 
the discrepancy demonstrated in the perturbative regime 
will be exactly compensated by higher order terms in the 
perturbative expansion. We will present a new result in 
the end of this paper which can be of relevance for ad- 
dressing the discrepancy. 



u> is the field frequency, top the plasma frequency, Ap the 
plasma wavelength and kp the plasma wavevector. 




FIG. 1: Parallel corrugated surfaces, with L representing the 
mean separation distance, a\ and ai the corrugation ampli- 
tudes and b the lateral mismatch between the crests. The 
corrugation are the smallest length scales in the perturbative 
expansion used in the paper, they have been exaggerated for 
the sake of a better visualization. 

In the following, we will suppose that the corrugation 
amplitudes are smaller than the other length scales 



01, a 2 < Ac, Ap, L 



(3) 



Within the range of validity of the PFA [19( , the Casimir 
energy in the presence of corrugations is simply obtained 
by adding over the mirrors' surfaces the contributions 
calculated with the local distances C 



II. GENERAL OUTLINE AND ASSUMPTIONS 

We first consider two parallel plane mirrors, Ml and 
M2, with corrugated surfaces (the case of a plane and a 
sphere with corrugations will be studied later on). The 
profiles of the mirrors Ml and M2 are defined by two 
functions h\{x, y) and hi(x, y) describing the local height 
with respect to mean planes Z\ = and zi = L. The 
mean planes are defined so that h\ and hi have null spa- 
tial averages, L representing the mean distance between 
the two surfaces; hi and hi are both counted as positive 
when they correspond to separation decreases. 

We assume that uniaxial sinusoidal corrugations are 
imprinted on both plates with the same period Ac and 
along the same direction, but with a spatial mismatch b 
between the corrugations crests (see Fig. [1} 



hi = a± cos(/cca ; ) > h% = a% cos (kc{x — b)) 
2tt 



kc = t— 
Ar 



(1) 



Ac is the corrugation wavelength and kc the correspond- 
ing wave vector. When a specific model of the optical 
response of the metallic mirrors will be needed, we will 
take the plasma model with the dielectric function 



e(u;) 



1 



2tt 
Ap" 
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1 Epp (£(r)) 
£pfa = / d r — 



C = L-hi-h 2 (4) 



Here Epp /A is the Casimir energy per unit area calcu- 
lated for plane and parallel plates, at the local separation 
distance £(r). Using the condition ([3]) and expanding (j4|) 
up to second order, we find the lowest-order correction 
to energy due to the presence of corrugations 
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h hi)' (5) 
aidi cos(kcb) 



As the energy corrections proportional to a 2 and a| do 
not depend on the lateral mismatch 6, they do not con- 
tribute to the lateral force 



TTilat 
^PFA 



dSSppA 1 d 2 Epp 



db 



2 dL 2 



kcaia 2 sin(fc c 6) (6) 



In the following, the expression of the force will be 
extended beyond the regime of validity of the PFA b 
using the general scattering approach presented in 
Performing a perturbative expansion up to the second 
order in the corrugation amplitudes, we will be able to 
evaluate the lowest-order energy correction due to the 
corrugations, proportional as to aiai cos(fccfr), but 
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with a different coefficient of proportionality in front of 
this quantity. We will then find a lateral force 



F lat = r PP ai a 2 sm(k C b) 



(7) 



with Tpp a function of the three length scales L, Ac and 
Ap. The main aim of this paper is to obtain the explicit 
expression of Tpp and to discuss various limiting cases. 

First, we will discuss in some detail the PFA limit, 
that is the case where (0 can be reduced to (J5]). This 
ap pro ximation was used for comparison with experiments 
in [23| and it has the advantage of being easily extended 
to higher orders in the corrugation amplitudes. However, 
as emphasized in [lj| [2(| for the case of roughness 
and then in [28| for the case of the lateral force, it can 
be valid only when the corrugation wavelength is larger 
than the other length scales (Ac S> L,Xp). A second 
interesting limit corresponds to perfect reflection of the 
mirrors. Then, the lateral force can be calculated in the 
path-integral theory developed in Ref. [24j, and higher 
order terms may also be calculated [25j]. We will prove 
below that the path-integral theory provides us with an 
independent test of our formalism since we recover its 
results in the limit Ap <C L, Ac- 

A third limiting case studied in the paper corresponds 
to the case of rugged corrugations Ac <C L, Ap. This case 
corresponds to evaluations far beyond the PFA regime 
and is particularly interesting as it constitutes a non 
trivial interplay between geometry and the Casimir ef- 
fect [ill [T3| • It is also of great interest for applications 
to the configurations with nanometric corrugations. We 
will derive analytical expressions in this limit and discuss 
the large deviation from PFA thus obtained. 

In the present paper, we will restrict our attention on 
the specific geometry where the corrugations of the two 
plates are aligned. The scattering approach also allows 
one to consider the case where the corrugations are not 
aligned so that the Casimir energy depends on the angle 
between the two directions, which results in a torque [3 11 ] - 
The measurement of the Casimir torque with torsion bal- 
ance techniques could be an alternative manner of testing 
the nontrivial geometry dependence of the Casimir ef- 
fect. Other surface profiles than sinusoidal corrugations 
can also be considered, as long as the amplitudes remain 
smaller than the other length scales. 



III. THE LATERAL CASIMIR FORCE 
BETWEEN CORRUGATED PLATES 

In this section, we derive a general expression for the 
Casimir energy between corrugated plates up to second 
order in the corrugation amplitudes. We start from a 
general expression for the Casimir energy valid in the case 
of arbitrary nonspecular scattering 0. This expression 
contains second-order correction terms proportional to 



a\ and a 2 , which were studied in detail in [2l|, but do 
not contribute to the lateral force. Here we focus on 
the second-order correction terms proportional to a\a 2l 
which are responsible for the lateral force. 

The electromagnetic fields are developed over Fourier 
components labeled by the two-dimensional wave vector 
k parallel to the r = (cc, y) plane and the polarization 
p (transverse electric, TE, or transverse magnetic, TM). 
Their scattering upon the non plane mirrors is then de- 
scribed in terms of non specular reflection operators cou- 
pling different wave vectors and polarizations (more de- 
tails in [H|). As a consequence of stationarity, scattering 
preserves the frequency lo. There exist two relevant non 
specular reflection operators, 1Z\ and IZ2 which describe 
respectively the intracavity fields reflected by the mirrors 
Ml and M2 as functions of the intracavity fields imping- 
ing these mirrors. 

The Casimir energy between the two non plane mirrors 
is then written as an integral over imaginary frequencies 
£ = -iu) [| 



£ 
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Trln(l-fti(^)e-^ft 2 (^) 



-KL\ 



(8) 



/C is a diagonal operator in the basis of plane waves (k, p) 
with the diagonal elements given by k = v/k 2 +£ 2 - The 
trace in Eq. ([5]) is a sum over the plane waves defined as 
/ pL) 2 I n the simplest configuration with two plane 
parallel plates, the reflection operator becomes diagonal 
with the diagonal elements given by the specular reflec- 
tion coefficients r-j ;p (k, We thus recover from (JSj> the 
known expression of the Casimir energy between plane 
parallel plates described by reflection amplitudes [33| 



HA 
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, 2nJ (2tt) 
d p (k) ee 1 - ri ;p (k)r 2 ;p(k) e 



V 

-2kL 



(9) 



The area A of plates has been introduced as a substitute 
for the quantity (2tt) 2 6^(0). 

The case of interest in this paper corresponds to cor- 
rugated plates. The reflection operators IZj thus contain 

zeroth-order contributions, denoted Hj and correspond- 
ing to plane surfaces, and non specular contributions in- 
duced by the corrugations 



(0) 



(10) 



Non specular operators 51Zj have been expanded in pow- 
ers of the Fourier transforms Hj (k) of the surface profiles 
hj (r) . It is now straightforward to write a perturbative 
expansion of the Casimir energy J5J) in powers of the cor- 
rugation amplitudes. The modification of the Casimir 
energy due to corrugations is read as 
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I? is a diagonal matrix with elements d p (k). We then collect the different orders in Hj. The first-order terms, which 
represent the change of Casimir energy due to the mean displacement of the mirrors, vanish thanks to the assumption 
(hj) = 0. The lowest-order corrections appearing at second order may be gathered in two categories, with square 
terms proportional to H\ or iJf (first line below) and cross terms proportional to H1H2 (second line), 



<5£pp = 



-KL 



(12) 



[7 + 1] denotes a sum modulo 2, and thus substitutes each 
mirror by the other. The second line has been simplified 
by noting that d~ 1 ri- p r2- p e~ 2nL = d^ 1 — 1, and then us- 
ing the invariance of the trace under cyclic permutations. 

The square terms (first line in Eq [T2"|) reproduce the 
roughness correction to the normal Casimir force, in full 
agreement with results reported in [2(| [2l[ . As they do 
not contribute to the lateral force, we disregard these 
terms in the sequel of the paper. In contrast, we focus 
our attention on the cross terms (second line in EqfT2"l) 
which generate the lateral force evaluated up to the sec- 
ond order in the corrugation amplitudes. We may in- 
cidentally note that the evaluation of the cross term 
is somewhat simpler than that of square terms, since 
the former is completely determined by the first-order 
non specular reflection operators 5TZj . The matrix ele- 
ments of these operators are simply proportional to the 
Fourier components of the surface profiles, with the form 
i2^,(k,k')flj(k-k'). The coefficients R^, depend on 
the optical properties of the mirror j and will be calcu- 
lated below for metallic mirrors described by the plasma 
model. 

Assuming that the two mirrors are made of the same 
medium, we write the cross contribution in (|12[) as 

r°° df f d 2 k' 

_ v e-(" , +")^(k^,k)flg, ) ,(k,kQ 

k '' k ^ dP(k)dP'(k') 

p' ,p 

For uni- axial sinusoidal corrugations (Q}, this second- 
order cross contribution to energy is read as 

^ P oss = I Q{kc) aia 2 cos(k c b) (14) 

Since Q(k) is negative in (|13|) . the energy is minimized 
when the crests are facing each other (6 = or a mul- 



tiple of Ac)- Differentiating with respect to the lateral 
mismatch 6, we obtain the expected expression ([7]) of the 
lateral force with the function Tpp defined by 

r PP = -g(k c )kc (15) 

In the next two sections, we check out that the known 
cases of PFA and perfect reflection are recovered as ap- 
propriate limits of our more general scattering expression. 
We then present more explicit results for the specific case 
of metallic mirrors. 



IV. THE PROXIMITY FORCE 
APPROXIMATION 

The Proximity Force Approximation (PFA) has been 
written as Eq. ([5]) above. In the context of our calcu- 
lations, we expect this expression to be recovered in the 
limit of very smooth surfaces Ac — > 00, that is precisely in 
the limit kcL, fee Ap <C 1. The question is thus whether 
or not the response function Q satisfies the condition 

The fact that this question has a positive answer can be 
proven in the context of our scattering formalism (2l| . 

To this purpose, we have to study the specular limit 
of the non-specular scattering formalism. In this specu- 
lar limit, the generalized reflection coefficients show the 
following behavior 

Jain R pp ] (k, k') = 2n r p 8 pp , (17) 

This can be checked out on the explicit expressions given 
below for the plasma model and is also true regardless of 
the model considered for the mirrors. As a matter of fact, 
Rj :pp i (k, k) would give the correction of the Casimir en- 
ergy for a mean displacement of the mirrors (hj) (should 
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the latter not be supposed to vanish). Hence it can be 
deduced from a global translation of the surface by a 
quantity (hj). For real frequencies, this amounts to an 
additional round-trip phase factor which is finally read 
as expression (1171) when going to imaginary frequencies. 

From Eq. (fT5|) , we now deduce the specular limit of the 
response function Q 



lim g(k) = -AH 

fc->0 



2tt 



d 2 k 



(2tt) 2 ^ 



(18) 



This expression turns out to fit condition (|16p when using 
the expression (J9j) of the Casimir energy Epp evaluated 
between plane plates. The property (|16p can be properly 
named as the "Proximity Force Theorem" with a simpli- 
fied expression indeed obtained in the secular limit of the 
more general non specular scattering theory. It holds in 
situations where specular scattering is sufficient to cal- 
culate the Casimir energy, that is also when the surfaces 
are approximately flat over distances of the order of L, 
since the main contributions to the Casimir effect come 
from wavelengths of the order of L. 

At this point, we want to emphasize that this dis- 
cussion does by no means imply that the limit Q16p 
can be used as an approximation for arbitrary values of 
fee- In particular, the lateral force computed within the 
PFA grows linearly with the corrugation wavelength kc 
whereas the scattering theory leads to a different behav- 
ior, eventually decreasing exponentially for large values 
of kc- In order to discuss deviation from PFA in an as 
clear as possible manner, we introduce the ratio between 
the force calculated in the scattering and PFA theories, 
in both cases at the lowest order oc a\a,2, 



P(kc) 



Tpp Q(k c ) 
S(0) 



pPFA 
1 pp 



(19) 



As a consequence of the discussions of the present section, 
this ratio goes to unity in the PFA limit 



lim p(k c ) 

kc— >0 



1 



(20) 



The variation of this ratio p for fcc / analyzed in the 
foregoing sections, measures the inaccuracy of the PFA. 



The effect of geometry on the Casimir effect between 
perfectly reflecting mirrors has been studied by various 
theoretical approaches (see [H, [TtJ for a review and 
[H, [35[ for recent results). The case of parallel plates 
with uni-axial corrugations has been studied in the per- 
turbative second-order approximation [24| as well as in 
the nonperturbative case [25|]. We restrict our attention 
here to the perturbative approximation, the results of 
which are reproduced exactly by our scattering formal- 
ism when the non specular reflection amplitudes known 
for perfect reflectors (see the Appendix C of Ref. [2l|) are 
plugged into the expression (fT5|) of the response function 
Q{k c ). 

According to the general method discussed in the pre- 
ceding section, we present the results in terms of the ratio 
(TT9"1) . Changing the integration variables to 7 = kL and 
7' — k! L, we obtain the following expression 



30 



TT^kr'L 



/ d7'e~ 7 ~ 7 

J\"f-k c L\ 

j[7 2 +7 ,2 -(fcc£) 2 ] 2 +7 2 7 /2 
(1 -e- 27 )(l-e- 27 ') 



(21) 



The functional dependence of p is very simple since it de- 
pends only on the dimcnsionless variable kcL = 2irL/\c 
which quantifies the smoothness of the surfaces on the 
scale determined by L. 

The PFA result p ~ 1 is recovered at the limit of small 
wavevectors kc —> (as predicted more generally by 
Eq. |20|) . In the opposite limit of large values for kcL, 
a completely different behavior is obtained with an ex- 
ponential decrease 



P : 



30 
^4 



(kcL) 4 
15 



(k c L) 2 + 3k c L + 3 



,-k C L 



(22) 



This expression is in full agreement with the results ob- 
tained in [24[ and this is also the case for the numerical 
evaluation of (|2ip represented on Fig. O We note that 
the PFA result is reproduced with some accuracy only 
for very small values of the dimensionless variable kcL. 

We now turn to the study of metallic mirrors described 
by the plasma model where there is a third length scale, 
the plasma wavelength Ap. We recall the expressions of 
the specular Fresnel reflection amplitudes 



V. THE LIMIT OF PERFECT REFLECTION 
AND THE PLASMA MODEL 

From now on, we take the plasma model ^ to de- 
scribe metallic mirrors. The ideal case of perfect reflec- 
tion is expected to be reproduced at the limit of very 
small plasma wavelength Ap <C L,\q. We first recall 
the results already known for ideal perfect reflectors and 
which only depend on the parameters L and Ac ■ We then 
apply the general formalism derived above to the plasma 
model and prove that the ideal case of perfect reflection 
is indeed reproduced when Ap is small enough. 



r TE (k,0 = 



k - Kt 

K + K t 



en — n* , 
r TM (k,C) = — (23) 



€K + K t 



e 7J = 



e = 1 



Kt denotes the imaginary part of the z component of the 
wavevector inside the metallic medium. We also intro- 
duce the following shorthand notations 



0=- , ft 

K 

h p (k,o 



k 

Kt 



P± 



K ± K t 
l±PPt 



1 -r p (k,£) 2 e- 2KL 



(24) 
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k c L 

FIG. 2: Variation of p versus the dimensionless variable kcL 
for the ideal case of perfect reflection. 

We then use the nonspecular first-order reflection coeffi- 
cients RpJ, (k, k') computed in [2l| with the help of the 
perturbation approach of (36|, based on the extinction 
theorem and the Rayleigh hypothesis [13, 0, [3§| . When 
inserted into Eq. (fT3|) . these expressions lead to 

&k,k' = ^2 HeV-'e' ( 25 ) 
e,e'=+,- 

x (/ lT E(k)/i T E(k')C 2 (1 + e/3/30 (1 + e'/?'#) 

-hTE(k.)hTM{k')S 2 (l + em) 
-/i T M(k)/iTE(k')5 2 (l + e'/3'/90 

+ h TM (k)h TM (k') (C + efifi) (C + e'/3'A)) 

The dependence on £ has been omitted, C = k • k' /(kk 1 ) 
and 5 = vl — C 2 represent the cosine and sine of the 
angle between the transverse vectors k and k'. 

A number of interesting properties can be checked out 
analytically on the expression obtained by plugging (|25p 
into (fT5)) . In particular, the PFA and perfectly-reflecting 
limits can again be recovered from this expression, by 
taking respectively kc <C kp, 1/L and fcp ^> fee, 1/L. Wc 
chose here to illustrate the same properties by discussing 
numerically integrated results. We plot p as a function 
of kcL, with fixed values of the second dimensionless 
quantity kpL. The values chosen for Fig. [3] ( fcp L = 1, 2.5, 
5 and 10) correspond respectively to L = 21.8, 54.5, 109, 
218 nm when taking Ap = 137 nm, which corresponds to 
gold-covered plates 41]. In all cases, p is smaller than 
unity and decreases when fcc increases. As already seen 
for perfect mirrors, the accuracy of the PFA is poorer 
and poorer for shorter corrugation wavelengths. At large 
values of fcpL, the curve obtained for the plasma model 
tends towards the curve calculated for perfect mirrors, 
as expected (compare the solid curve on Fig. [3] with that 
on Fig. [2]). Otherwise, the result obtained for the plasma 
model differs from that for perfect mirrors. 

Let us give a few numbers here, with parameters 



0.0' ' 

1 2 3 4 5 

k c L 

FIG. 3: Variation of p versus the dimensionless variable 
kcL for metallic mirrors described by the plasma model, for 
kpL =1 (dashed line), 2.5 (dotted line), 5 (dashed-dotted line) 
and 10 (solid line) [colors online with respectively green, blue, 
red and black lines]. 



L = 220 nm, Ac = 1.2 pm and Ap = 137 nm chosen to 
be close to the experimental figures of (23[. As stated 
in the preceding paragraph, this corresponds to a large 
value kpL = 2irL/\p ~ 10, so that the value of p cal- 
culated for the plasma model approaches that obtained 
for perfect mirrors. The precise values are p = 0.814 for 
the plasma model and p = 0.819 for perfect mirrors. The 
important point to be noticed here is that, contrarily to 
some claims [29j, these values lie far from the PFA ex- 
pectation (which is simply p = 1). The discrepancy is 
larger in the experiment which employs a plane-sphere 
setup, as discussed in more detail below. 

To make this point completely clear, it is true that, 
for perfect mirrors, p is determined exclusively by the 
quantity k^L (see discussions above). This is still ap- 
proximately true for the plasma model as soon as the 
second dimensional quantity fcpL is large. Nevertheless, 
this property cannot be considered as universal since the 
function Tpp is generally a function of the three variables 
L, fcc and fcp, and the dimensionless function p a func- 
tion of the two dimensionless variables fccL and fcpL. To 
give an example, L — 55 nm and Ac = 300 nm lead to 
the same value of fcc^ as in the preceding paragraph. 
The values p = 0.838 found for the plasma model thus 
significantly differs from p = 0.819 for perfect mirrors. 
Even more striking illustrations of the dependence of p 
on fcpi will be given in the next section where we discuss 
the regime far beyond the PFA. 



VI. RUGGED CORRUGATION LIMIT 

We now consider the limit of short corrugation wave- 
lengths Ac <C L, Ap which, as already noticed, corre- 
sponds to a non trivial interplay between geometry and 
the Casimir effect [l|| [ItJ while being of great interest 
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for the problem of nanostructured surfaces. 

We have seen that p decreases exponentially with 
fcci > 1, a property which can be given a simple ex- 
planation [2l|, |40J. The cross contribution to corruga- 
tion energy fT^]) , evaluated at the lowest order, origi- 
nates from intracavity propagation loops containing one 
non-specular reflection at each mirror. The first non- 
specular reflection changes the field momentum from k' 
to k, whereas the second changes the momentum back 
to its initial value (see Eq. |2"5]) . The exponential fac- 
tors exp(— kL) and exp(— k'L) associated with intracav- 
ity propagation (for imaginary frequencies) thus lead to 
the following behavior in the rugged corrugation limit 

p~[3{k c L) 7 / 2 e- kcL , Ac«L,A P (26) 




k c L 



(3 is independent of kcL while remaining a function of 
kpL. 

Equation (|26p corresponds to a behavior very different 
from that obtained for perfect mirrors (see Eq. |2"2"]) . This 
entails that p takes very different values at the same kcL 
but different kpL, with 



p(k c L > l,k c L > k P L) 
p(k c L > l,k c L < k P L) 



(k C L) 



-1/2 



< 1 



(27) 



In other words, the perfectly-reflecting model grossly 
overestimates the lateral Casimir force obtained in the 
limit of rugged corrugations. Note that this occurs for a 
large separation distance L though the latter condition is 
often carelessly associated with perfect reflection. This 
is a serious limitation of the model of perfect reflection 
which, as already emphasized, can only be trusted when 
Ap is much smaller than Ac and L. 

To propose a better visualization of this feature, we 
plot on Fig. 0] the quantity a = pexp(kcL) for a large 
distance L = 1 pm (kpL = 46 with Ap =137nm), as a 
function of k^L. The result shown as the solid curve 
differs almost everywhere from the result obtained for 
perfect mirrors a = (kcL) 4 , shown as the dashed line. 
As kcL increases past kpL, a stays below the result for 
perfect reflectors. 

The exponential fall-off of p at large kcL has to be 
contrasted with the linear increase obtained in previous 
papers for the similar function associated with roughness 
[2l|. Such a linear growth was resulting from the con- 
tribution of second-order non specular reflections which 
do not contribute to the lateral force. It follows that, 
whereas the PFA was underestimating the roughness cor- 
rection, it overestimates the lateral Casimir force. 



VII. COMPARISON WITH EXPERIMENTS: 
THE PLANE-SPHERE CASE 

We now consider the lateral Casimir force in the plane- 
sphere (PS) configuration, which corresponds to the ex- 
periments [231 ] . We derive the PS result from the plane- 
plane (PP) one by using the PFA for treating the sphere 



FIG. 4: Variation of a versus kcL (solid line); the high-fc limit 
for perfect reflectors is shown as the dashed line [colors online 
with respectively black and purple lines] . 



curvature effect. The validity conditions required here 
are much more easily met than those needed for the cor- 
rugation effect. Indeed, experiments use spheres with 
radius R of hundreds of microns, of the order of thou- 
sand times larger than the distance L. The condition 
RL ^> is also needed in order to treat curvature and 
corrugation effects without taking any complicated inter- 
dependence into account, and it is met with the experi- 
mental figures [23l ]. 

We then use the PFA to deal with the sphere curvature 
effect, while accounting for deviations from the PFA for 
the corrugation effect. We obtain the cross energy cor- 
rection S£pg ss in the PS case from the one 8£pp ss calcu- 
lated for the PP case by integrating over the distribution 
of distances generated by the sphere curvature 



<5£p r s oss 



CM) 



* rc cross 

2nRdL' 



(L',b) 



(28) 



b is the lateral mismatch of the parallel corrugations on 
the two plates and L the distance of closest approach. 
When differentiating with respect to b, one obtains the 
lateral Casimir force in the PS configuration as 



8 f°° 
F P a i = --S£™f(L,b) = J 2nRdL 



(29) 



This is then read exactly as expression ([7]), with a func- 
tion Tps now written as (using Eq[l5| 



/>oo 

r PS = J nkcRdL' g L ,{k c 



(30) 



Should the PFA be applied to the corrugation effect at all 
distances L' > L involved in this integral (|30|) . GL'(kc) 
would be replaced by Gl'(0) (see Eq. [16]), leading to 



f pfa , n f PP 
1 pg = tt/cc-K— 



d£i 



pp 



"pp 



dL 



(31) 
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The accuracy of this approximation is well represented 
by the following ratio, defined by analogy with (fT9"|) , 



pps(fcc) 



i PS 
pPFA 
1 PS 



Fpp(i) 



(32) 



dL' 



-dF PP (L') 
dL> 



PL>{kc) 



It is an immediate consequence of this expression that 
the PFA accuracy will be worse in the PS configuration 
than in the PP one (for the same L) since it is determined 
by the PP case with distances L' larger than L. 




FIG. 5: Lateral force coefficient Tps (force divided by 0,10,2) 
for the plane-sphere geometry, as a function of kcL (L and 
Ap chosen to fit [2^]); the solid line is the result of scattering 
theory while the dotted line corresponds to the PFA expres- 
sion; the value kcL met in experiments is shown as the verti- 
cal dashed line [colors online with respectively black, red and 
blue lines]. 

In order to stay closer to the experimental figures, we 
have chosen here to illustrate these results by plotting 
on Fig. [5] Tps as a function of kcL. The parameters 
L = 220nm and Ap =I37nm are chosen to fit numbers 
of [23[ . The solid line representing the result of scat- 
tering theory has to be compared with the dotted line 
associated with the PFA. The experimental corrugation 
wavelength kc — 5.2 pm^ 1 is indicated by the vertical 
dashed line. It is clear that the value thus attained by 
k^L =1.14 is large enough to produce a significant inac- 
curacy of the PFA. The precise numbers for this value 
are r PS = 421 pN/^m 2 and rp| A = 585pN/^m 2 respec- 
tively for scattering theory and PFA. In other words, the 
scattering theory result is smaller by a factor pps — 72% 
than the PFA result with the same choice of parameters. 
As expected, the inaccuracy of the PFA is significantly 
worse than it was in the PP case where p = 81% was 
found (with the same choice of parameters). 

Let us now discuss Fig. [5] in a more general manner. 
As previously, the PFA result (f3Tj) grows linearly with 
kc and this is also true for the scattering theory result at 
small kcL. When kcL increases, the scattering theory 
result begins to deviate from PFA and eventually decays 



exponentially at large values of kcL. Tp§ thus shows a 
peak value, found to lie at kcL — 2.08 (that is Ac = 
665 nm when L = 220 nm). For the PP setup, the peak 
position was found at k<^L = 2.6 and the difference can 
again be explained from the fact that the PS force is an 
average of the PP result over L' > L. 

We have also plotted on Fig. [6] Tps as a function of 
L with Ac = 1.2 /xm and Ap = 137 nm fixed at their 
experimental values. The lateral Casimir force is found 
to decrease with the distance L. The scattering theory 
result (solid line) corresponds to an exponential fall-off 
for L > Ac whereas the PFA result (dot-dashed line) 
decays only as a power law, L~ 3 in the plasmon regime 
(L <C Ap, dotted line) and L~ 4 in the perfectly-reflecting 
regime (L 3> Ap, dashed line). The scattering theory and 
PFA results agree at short separation distances L < Ac. 
For very short distances L <C Ap < Ac furthermore, the 
force can be deduced from the plasmon approximation 
[13, 511]. The crossover between the L~ 3 behavior in this 
range and the exponential decay at long distances occurs 
in the region around a few hundred nanometers, which is 
magnified in the inset of Fig. [6] Note that this is also the 
range tested experimentally, which we will discuss now. 




10 100 

L (nm) 



1000 



FIG. 6: Lateral force amplitude coefficient Fps (force divided 
by a\a2) in the plane-sphere geometry, as a function of L, 
with Ac and Ap chosen to fit [23]; the curves correspond to 
scattering theory (solid line), PFA (dotted-dashed line), PFA 
combined with perfect reflection (dashed line) and plasmon 
(dotted line) approximations [colors online with respectively 
black, red, blue and green lines]. 

We come now to a discussion of the number ppg ~ 72%, 
which measures the inaccuracy of PFA, and points at 
a potential concern for the theory-experiment compari- 
son. We recall that the experiments were performed with 
corrugation amplitudes smaller, but not much smaller 
than the other length scales ([23| report a\ = 59 nm, 
02 = 8nm, to be compared to Ap = 137 nm, L = 
220 nm, Ac = 1.2 pm). As already discussed in [28|], this 
point made a direct comparison between experiments and 
our perturbative scattering theory impossible, and thus 
pushed us to try to discard the higher orders contribu- 
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tion to the PFA calculation of [23[ in order to perform 
an indirect comparison. As higher order corrections have 
not yet been estimated within scattering theory, we will 
not pursue this line of reasoning further. 

Here we want to emphasize the comparison between 
results obtained within PFA and beyond PFA for the fig- 
ures of experimental interest. When restricting the atten- 
tion to calculations up to the second order, this compar- 
ison is precisely characterized by the number pps — 72% 
calculated for the parameters favored in |29j ] . In our opin- 
ion, this discrepancy (~ 28%) does not lie so far from the 
margins of experimental uncertainty (0.32±0.077pN in 
[1^]), which correspond to a relative accuracy of ±24%. 
In contrast to this opinion, the concern was made more 
acute by a Comment [29[ which fabricated a larger dis- 
crepancy by comparing two numbers which are not to be 
compared (and which we did not compare), namely the 
perturbative result beyond the PFA and the non pertur- 
bative result within the PFA. We have already explained 
our point of view in a Reply [30l |. and now present new 
discussions of the issue in the sequel of the section. 

An interesting way of addressing the issue is suggested 
by a close scrutiny of Fig. [6] If we fit the result of scat- 
tering theory to a power law within the interval of exper- 
imental interest, we obtain a law oc L~ AA in agreement 
with experiment [23| . This coincidence is also apparent 
in the fact that the scattering theory curve is roughly 
parallel in this distance interval to the law obtained by 
combining PFA and the model of perfect reflection (see 
the inset of Fig. [6]). This suggests that the difference 
between these two curves can be confused with a poor 
determination of the absolute value of the distance L 
between the mirrors. As far as theory is concerned, L 
is precisely defined as the mean separation distance be- 
tween the corrugated surfaces, such that the corrugation 
profiles have zero mean values (see Fig. [1} . Precise mea- 
surements of variations of L were reported in [23| , but the 
absolute determination of its value was by far more dif- 
ficult, in a geometrical configuration where not so small 
corrugations are facing each other on plane and spherical 
surfaces. 

In order to test the idea that the difference between 
scattering theory and PFA could be confused with an 
offset in the determination of the absolute distance L, we 
plot on Fig. [7] three curves in the range from 220 to 260 
nm. The solid and dotted lines represent the scattering 
theory and PFA results calculated with the value of L 
supposed to be ideally determined; the third, dashed- 
dotted, line describes the scattering theory result with an 
offset of 20 nm for the distance L. The important point 
is that the scattering theory with an offset can easily be 
confused with the PFA without it. More precisely, the 
difference between the two curves is certainly within the 
margins of experimental uncertainty, the magnitude of 
which is of the order of the initial difference between the 
scattering theory and PFA results. 

This means that the difference between scattering the- 
ory and PFA can have been confused with an offset of 20 
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l -I 

220 230 240 250 260 

L (nm) 

FIG. 7: Same conventions as in Fig. [5] with curves restricted 
to the range of distances from 220 to 260 nm; the solid and 
dashed-dotted lines represent the scattering theory and PFA 
results calculated for the value of L indicated by the abscissa; 
the dotted line is the scattering result computed with an offset 
of 20 nm for L [colors online with respectively black, red and 
blue lines]. 

nm in the distance measurement. Note that this value 
is well below the larger corrugation amplitude (59 nm) 
in the experiment. It is therefore consistent with the al- 
ready discussed experimental difficulty for localizing the 
position of the abstract reference planes at a scale smaller 
than the corrugation amplitudes. Note also that, if the 
absolute determination of the distance has been helped 
by a global fit of the experimental results to the theory, 
the offset should have been automatically produced by 
the fitting procedure. 

VIII. CONCLUSION 

We have studied the lateral Casimir force arising be- 
tween two corrugated metallic plates, using scattering 
theory in a perturbative expansion valid when the cor- 
rugation amplitudes are smaller than the other length 
scales L, Ac and Ap. We have shown that the Proximity 
Force Approximation (PFA) is recovered at the limit of 
smooth plates Ac ~> L, Ap and we have also obtained an 
expression for the lateral force in the opposite limit of 
rugged corrugations Ac <C L, Ap. We have reproduced 
the results known for perfect mirrors when Ap <C Ac, L 
and have also given expressions valid when this is not the 
case. 

As the perturbation conditions a\,a2 <C £,Ac,Ap 
are not met in the experiment [23| . the comparison 
of scattering theory with measurements of the lateral 
force will require more work. Progress on this question 
could be achieved by calculating higher order corrections 
for metallic mirrors beyond the PFA. These corrections 
would affect the numbers given in the present paper, but 
it is unlikely that they could compensate exactly the de- 
viation from PFA which has been demonstrated in the 
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perturbative theory. 

Progress could alternatively come from further experi- 
ments. It would be very useful to have experiments with 
small corrugation amplitudes attaining the domain of va- 
lidity of the perturbative theory. Other relevant improve- 
ments would be to increase the experimental accuracy, in 
order to be able to distinguish more easily between alter- 
native theories, and if possible to measure the absolute 
distance in a more reliable manner, in order to get rid of 
the offset confusion. Of course, this program raises seri- 
ous experimental challenges, given the minuteness of the 
lateral force effect. But the reward would be remarkable 
with a potential experimental access to a configuration 



where arises a nontrivial interplay between geometry and 
the Casimir effect, that is beyond the PFA. 
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